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Introduction 
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At short distances QCD vacuum mimics the Fock space ground state of pertur- 
bation theory: deep inelastic scattering experiments, jet production in high energy 
P> '. reactions, QCD sum rules provide empirical evidence for that. 

In fact colour is confined, and quarks and gluons never appear as free particles in 
asymptotic states: the ground state is very different from the perturbative vacuum. Its 
exact structure is not known. Many models have been attempted to give an approx- 
imate description of it: some of them are based on "mechanisms", i.e. they assume 
that the degrees of freedom relevant for large distance physics behave as other well 
^ . understood physical systems. 

The most attractive mechanism for colour confinement is Dual superconductivity 
of type II of QCD vacuum 1, 2 Dual means interchange of electric with magnetic with 
respect to ordinary superconductors. The idea is that the chromoelectric field in the 
region of space between a QQ pair is constrained by dual Meissner effect into Abrikosov 3 
flux tubes, with constant energy per unit length. The energy is then proportional to 
the distance 

E = oR (1) 

and this means confinement. 

Lattice is the ideal tool to study (at least numerically) large distance phenomena 
from first principles. There is indeed evidence from lattice simulations that: 

1) The string tension exists: large Wilson loops W(R,T) describing a pair of static 
quarks at a distance R propagating for a time T, obey the area law 4 

W(R,T) ~ exp(-aRT) (2) 

Since in general W(R,T) ~ exp(— V(R)T), the observed behaviour Eq.@ con- 
firms confinement, as defined by Eq.(^). 
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2) Chromoelectric flux tubes have been observed, joining Q Q pair propagating in 
Wilson loops 5, 6 . Their transverse size is ~ 0.5 fm. 

3) String like modes of these flux tubes have been detected 7 . 

4) Particles belonging to higher representations than quarks also experience a string 
tension at intermediate distances, which, for SU(2), depends on the colour spin 

J as 8 

a j = kJ(J + l) 

or 

^- = 3 jJ(J+l) (3) 

01/2 4 

Observations l)-3) support the idea of dual superconductivity. We will discuss in detail 
the implications of 4) in the following. 

The problem we adress in these lectures is: can dual supeconductivity of QCD 
vacuum be directly tested? 

The plan of the lectures is as follows. We will recall basic superconductivity and its 
order parameter, in order to clarify what we are looking for. We will then define dual 
superconductivity and its disorder parameter. We will construct the disorder parameter 
for U(l) pure gauge theory. 

We will then check the construction with the X — Y 3d model (liquid He4). 

We will then revisit the abelian projection, which reduces the problem of dual 
superconductivity in QCD to a U{\) problem. The Heisenberg ferromagnet will prove 
a useful laboratory to check this procedure. We will finally show how dual supercon- 
ductivity can be directly detected in 577(2) and SU (3) gauge theories. 

A discussion of the results and of their physical consequences is contained in the 
final section. 



Basic superconductivity: the order parameter 9 . 

A relativistic version of a superconductor is the abelian Higgs model 

C = -\F, V F, U + (D^(D^) - V{$>) (4) 

is the electromagnetic field strength, is the covariant derivative 

D lt $ = (d lt -iqA li )$ (5) 
and V($) the potential of the scalar field 

V(*) = \(&*-v?) 2 (6) 

If /i 2 is positive the field $ has a nonzero vacuum expectation value. Since $ is a 
charged field this is nothing but a spontaneous breaking of the £7(1) symmetry related 
to charge conservation. The ground state is a superposition of states with different 
electric charges, a phenomenon which is usually called "condensation" of charges. 
A convenient parametrization of $ is 

$ = pe^ p = p t > o (7) 
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Under gauge transformations 



Ay — ► — d^a + a (8) 

The covariant derivative (|5|) reads in this notation 

= e Wc > [d, - iq(A, - d,6)} p (9) 

The quantity A tl = A^ — d^Q is gauge invariant. Moreover 

Fti,u = dp, A v - d v A^ = d»A v - d u A^ (10) 

The equation of motion reads, neglecting loop corrections, (or looking at C as an 
effective lagrangean) 

+ ^-A u = m = v / 2g($) (11) 

In the gauge A = a static configuration has d$A = 0, <9o$ = so that Ei = F 0i = 0. 
Eq . (|TT1) implies that 

V AH + —A = (12) 

The term rh 2 /2A in Eq. fll^) is a consequence of spontaneous symmetry breaking and 
is an electric stationary current (London current). A persistent current with E — 0, 
means p = since pj = E and hence superconductivity. 
The curl of Eq. (|12"D , reads 

o - m 2 - 

V 2 H H = (13) 



The magnetic field has a finite penetration depth 1/m, and this is nothing but Meissner 
effect. The key parameter is ($), which is the order parameter for superconductivity: 
it signals spontaneous breaking of charge conservation. 

Besides = 1/m there is another parameter with dimension of a length, A$ = 

If > V2X& the superconductor is called Ilnd kind otherwise it is 1st kind. 

For a superconductor of first kind,there is Meissner effect for external magnetic 
field H < H c (critical field), for H > H c the field penetrates the bulk and supercon- 
ductivity is destroyed. For second kind instead a penetration by Abrikosov flux tubes 
of transverse size 1/m is energetically favoured. Flux tubes repel each other. By in- 
creasing the external field the number of flux tubes increases. When they touch each 
other the field penetrates the bulk and superconductivity is destroyed. 

In a dual superconductor the role of electric and magnetic field is interchanged. 
The U(l) symmetry related to magnetic charge conservation is spontaneously broken, 
i.e. monopoles condense in the vacuum. An order parameter for dual superconductivity 
will then be the vacuum expectation value of a field carrying non zero magnetic charge. 

Monopoles. 

The equations of motion for the electromagnetic field in the presence of an electric 
current j M and of a magnetic current jjf are 

d^F^ = j v 

d"F; u = j« (14) 
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If both jp and jff are zero (no charges, no monopoles) photons are free, and the 
equations of motion are invariant under the transformation 

F^ u -> cos 9 F^ + sin 9 F* u 

F;„^cos9F; v -syn9F fll/ (15) 

for any 9. In particular if 9 = tt/2 Eq.'s(15) give E — > H, H — > — E which is known 
as duality transformation. In nature jjf = 0. The general solution of Eq.(|T3D is then 
written in terms of vector potential 

= d^A u - d v A ll (16) 

and d^F* = is identically satisfied (Bianchi identities). 

If a monopole exists Bianchi identities are violated. However they can be preserved, 
and with them the description in terms of A^, by considering the monopole as the end 
point of a thin solenoid (Dirac string) connecting it to infinity: the flux of the Coulomb 
like magnetic field, H = jf^g, &(H) = M is conveyed to infinity by the string 10 . 

The string is invisible if the parallel transport of any electric charge around it is 
trivial: 

Q iq§ Adx _ e iq$(H) _ ^jrn _ ^ 

or 

qM = 27m (17) 

Eq.flTT|) is known as Dirac quantization condition, and constrains the U(l) group to be 
compact. If one insists to describe the system in terms of A^ monopoles are non local 
objects with non trivial topology. One could introduce dual vector potential B^, such 
that F* = d^B u — d u B^. Dual Bianchi identities would read d^F^ = 0, monopoles 
would be pointlike but electric charges could only exist if dual strings were attached to 
them. 

There is another acceptation of duality, which originates from statistical mechanics. 
The prototype example is the 2d Ising model. The model is defined on a square lattice, 
by associating to each site i a field a(i) wich can assume 2 values, say ±1. The action 
can be written 

£ = -/?5>(iM;) (18) 

the sum running on nearest neighbours. The partition function K[f3] is known exactly 
in the thermodinamical limit. At high j3 (low temperatures) the system is magnetized 
(cr(i)) 7^ 0; at low (3 it is disordered. A dual description can be given of the same 
system, by associating to each link (dual lattice site) a variable a* with value —1 if the 
values of a in the sites connected by the link are the same, +1 otherwise. 

It can be rigorously proven that the partition function in terms of the new variables 
has the same form as the original one K[/3] 



K*[(3] — K[(3*\ (19) 



with the only change 



P^IT ,r = Iarcsinh(-y (20) 

A relation like Eq . (|19"D is called a duality relation. It maps high temperature (strong 
couplings) regimes of K* to the low temperature (weak coupling) of K 11 . 
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Similar relations have been recently discovered in SUSY QCD with N = 2, and 
more generally in models of string theory 12 . 





t 




X 



Fig.l A kink in Ising model. 



If we look at the Ising model as the euclidean version of some 1+1 dimensional field 
theory, a configuration at fixed t, with a(i) — — 1, i < i a(i) — +1, i > i will appear 
on the dual lattice as a single spin up. This configuration has topology, it is a kink. The 
exitations of the dual lattice are kinks. At low temperature the system is magnetized 
(a) ^ and very few kinks are present. At high temperature (a) = 0, but, by duality 
relation, (a*) ^ 0. (a*) is called a disorder parameter, as opposite to (a) which is the 
order parameter. The relation 



can be proven in the thermodynamical limit, (a*) signals the condensation of kinks. 

We shall next adress the study of monopole condensation in U(l) compact gauge 
theory. We shall define a disorder parameter for this system, which describes the 
condensation of monopoles in the vacuum at high temperature (low f3). The parameter 
will be the v.e.v. of an operator with non zero monopole charge, and will thus signal 
dual superconductivity. 

The same construction will then be used for non abelian gauge theories after 
abelian projection. 

Monopole condensation in compact U(l): a disorder parameter 13 . 

Like any other gauge theory, compact U(l) is defined in terms of parallel transport 
along the links joining nearest neighbours on the lattice 



a being the lattice spacing. 

In the following we shall denote eaA^n) as O^n). The action is written in terms 
of the parallel transport ILj^ around the elementary square of the lattice in the plane 




(21) 



Uf,(n) = exp(ieaA M (n)) 



(22) 




(23) 



(24) 



As a 








(25) 



4 



5 



and Eg. (p^D describes photons if /3 = 1/e 2 . 

The generating functional of the theory (partition function) is 



m = [ n 



d9 pip) 



2tt 



exp(-S) (26) 



The theory is compact, since S depends on the cos of the angular variables, and is 
invariant under change of variables 

with arbitrary /^(n). A special case of Eq.(^) are gauge transformations. 

A critical /3, (3 C ~ 1.011 exists such that for (3 > (3 C the theory describes free 
photons. For f3 < j3 c electric charge is confined: Wilson loop obey area law 14 Eq.(Q) 
and flux tubes are observed 15 . 

A variant of the theory is provided by the Villain action 



exp(-S) = J2 ex P 



■f J2 \e^ v -2^m\ 



(2? 



For this variant, condensation of monopoles has rigorously been proven as a mechanism 
of confinement 16 . Recently the proof has been extended to more general forms of the 
action, including Wilson action Eq. (|23|) 17 . Monopoles are identified and counted by the 
following procedure 18 . 

Since by construction tc < 0^{n) < it, it follows from Eq . (|23|) that 

-4vr < 9^{n) < 4tt 

Qpv can be redefined modulo an integer multiple of 2ir, 

V = 0»v + 27rn^ - 7r < 9^ < tt (29) 

and the monopole current as 

Pu = ^nvpa^un pa = d u F* fiU (30) 

The total number of monopoles is 

A^EPoV) (31) 

n 

N M is large in the confined phase, and drops to zero in the deconfined phase. It has 
sometime been identified with the disorder parameter for monopole condensation. Of 
course N M commutes with the monopole charge, and therefore cannot signal by any 
means spontaneous breaking of magnetic U(l). 

The disorder parameter. 

The basic idea of the construction 11 ' 19 of a disorder parameter is the simple formula 
for translations 

e^lar) = \x + a) 

If we identify in our field theory 

x — »■ A(x, t) 

p — > E(x, t) — — i— =r^— - (32) 
6A(x,t) 
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then the operator 



fjt(y,t) = exp 



d 3 xE(x, t)b(x - y) 



(33) 



operating on field states in the Schrddinger representation will give 

n\A[x,t)) = \A(x,t) + -b(x - y)) 

— * 

i.e. it will add a monopole to any field configuration provided that b is the vector 
potential describing the field produced by the monopole 



-b(f) 



1 m r An 
e 2-7T r(r — r ■ n) 



(34) 



The gauge has been chosen to have the Dirac string in the direction n. 

\i is independent on the choice of the gauge for b if E obeys Gauss law. 
On a lattice 20 , after Wick rotation, and with the identification 



a?E l = -Imir + 0(a 4 

e 



~ - sin 6 0i 

e 



(35) 



Mif,no) = ex P 



(36) 



Here bi(n) is the discretized transcription of Eq.(|3"4"D and the (3 in front comes from the 
normalization 1/e in Eq.fl2"2"|) times the 1/e appearing in the monopole charge. 

A better definition (compactified) which shifts the angle 8 0i and not sin # i is 13 



fi(y,n Q ) =exp(3j2 [S(0 w (tf,n o ) + b\n - y)) - S(6 0i (n,n )) 



(37) 



which reduces to Eq.(^6j) at first order in b. 

In Eq.(|37]) by S we denote the density of action. The action can be any form, e.g. 
Eq.(|24l) or ( p8|) provided Eq.(^) is satisfied. 

If any number of monopoles or antimonopoles are created at time no, then b l {n — y) 
has to be replaced by the corresponding field configuration. To compute correlation 
functions of operators at different times, the rule is 



(KVh . . . n(y k , n° k )) = - J [] 



2tt 



exp (P(S + AS)) 



(3f 



where S + AS is obtained by replacing in the action the plaquettes II 
cos(0 Oi (n>°J) by 1 



n, n 



•> '"a J 



COS 



10/ 



(n,n° a ) + b\n-y a )). (1 < a < k). 



In particular we will study the correlator 



V{x ) = (/i(0,x o )/i(0,0)) 



(39) 



At large enough x 



V{x ) ^ Aexp{-M\x \) + (fi) 2 

\xo\— > °o 



(40) 



The last equality follows from cluster property, translation invariance and C invariance 
of the vacuum. Our aim will be to extract M and (/i) from numerical determinations 
of T>(xq). (fi) is the disorder parameter: a non zero value of it in the thermodinamic 
limit signals dual superconductivity. M is the mass of the lightest excitation with the 
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quantum numbers of a monopole and is a lower limit to the mass of the effective Higgs 
field which produces superconductivity. 
As explained above 

V(x ) = ^-Z[S + AS] (41) 
where S + AS is obtained from S by the change 

e i0 {n, o) -> e i0 {n, o) + &*(n) (42) 

0*°(n, x ) -> io (n, x ) - b\n) (43) 
and 6 l (n) is defined by Eq.(|34"D. Since 

6 M (n, 0) = -6 l (n, 1) + 6\n, 0) + 0°(n + i, 0) - 0°(n, 0) 
the replacement Eq.(f4"2"|) amounts to the change 

0Xn, 1) -> ^(n, 1) - 6i(n) = 0*(n, 1) (44) 



The change Eq.(|44|) can be reabsorbed in a redefinition of variables of the Feynman 
integral defining Z[S + AS 1 ] which leaves the measure invariant [Eq.(27)]. As a conse- 
quence 

y (n, 1) -> ^'(n, 1) + AW - A^ 
meaning that a monopole is added at xq — 1. Moreover 

Ol (n,l) -> Oi (n,l) +6*(fi) 

Again this change can be reabsorbed by a change of variables 

0*(fi,2) -> 0*(n,2) -foi(n) 

after which 

0^'(n, 2) -> ij '(n, 2) + A l b j - A j b i 

and 

lO (n,2) -> iO (n,2) +6*(tJ) 

The construction can be repeated till n = x , when the addition of &j(n) to 0*°(n, n ) 
cancels with the term —bi(n) in Eq . fl43|) . The change Z[S + AS 1 ] amounts to add a 
monopole in the site n = y, propagating from time to time xq. 

Measuring V(xq), to extract M and (/i) is non trivial, due to large fluctuations. 
AS" is the change of the action on a spatial volume V: it fluctuates roughly as y/V, 
which means a fluctuation ~ exp(\/V) on (//). A way out of this difficulty is to 
measure, instead of V(x ), the quantity 

p(x ) = ^ In Vino) = (S) s - (S + AS) s+ ss (45) 

The last equality trivially follows from the definition of Z. (S)s means average of S 
computed by weighting with the action S. Since £>(xo)/3=o = 1 



V(x ) = exp N 



%(/WJ (46) 



As xo — > oo 

p ~ 2-^ln( yU )+Cexp(-Mx ) (47) 
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Fig. 2 shows a behaviour of p as a function of xq consistent with Eq. ((47|) . 
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Fig. 2 Monopole antimonopole correlation. (Lattice 8 3 x 16) 
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Fig. 3 shows p(oo) = 2^1n(/i) as a function of (3. A huge negative peak appears at the 



phase transition, which, according to the definitions Eq . (f45|) , Eq.fllq) reflects a sharp 
decrease of (fj). 
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Fie. 3 function of f3. The negative peak signals the phase transition. 

(Lattice 8 3 x 16) 



This can be appreciated from Fig.4 where a direct measurement of (/i) is shown, even 
with large errors. 
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Fig.4 (p) v.s. /3. Lattice 10 3 x 20. 
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Fig. 5 poo versus 1/L for (3 = 1.009. 

At (3 > (3 C the system describes free photons: p and p(oo) can be computed in pertur- 
bation thory by a gaussian integration. The numeric result is, for a lattice L 3 x 2L 

Poo = -10.1 • L + 9.542 (P > /3 C ) (48) 

Poo tends to — oo in the thermodynamical limit L — > oo, or (p) — >• 0. In fact (p) is an 
analytic function of /3 at finite volume, and cannot be exactly zero for (3 > j3 c , since it 
would be identically zero everywhere. Only as V — > oo Lee - Yang singularities develop 
and (p) = for f3 > (3 C ( |48| ) as a respectable order parameter. 
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For (3 < f3 c , poo tends to a finite value as V — > oo. This can be seen from fig. 5 but 
it is also a theorem proven in ref.(16) which generalizes the result of ref.(15) for Villain 
action. 

For f3 ~ f3 c there is a phase transition, which is weak first order or second order. In any 
case the correlation length will grow as (3 — > (3 C in a certain interval of (fl c — P) with 
some effective critical index v: 

S~{p c -P)- V (49) 
A finite size scaling analysis can be done as follow. In general, by dimensional reasons 

( L a\ 



As (3 — > (3 C , a/£ — > and p 



which implies in turn that 



n = n 

//(L/CO) or 
p = f(L 1 /»(P c 



(50) 



d/J W f(L ((3 C — 0)) 



(51) 



(52) 



or pL~ l l u is a universal function of L l l u (f3 c — (3). 

Data from lattices of different size will lay on the same universal curve only for 
appropriate values of f3 c and v. The best values can be then determined. We obtain 

P c = 1.01160(5) v = 0.29(2) (53) 

ji ~ (p c — P) 5 as P — > Pc, and 

P „ 5 
L Vv L^(p c -P) 

For 5 we obtain 5 = 1.1 ±0.2. 

Our result is then that (//) ^ for P < P c in the thermodinamic limit, i.e. that 
the system is a dual superconductor for P < P c . 

We have also measured the penetration depth of the electric field, i.e. the mass 
of the photon, m, by the method of ref.(18). m properly scales as P — > P c with index 
v, and the indication is that it is substantially smaller than M, fig. 6, or that the 
superconductor is type II. 
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Fig. 6 Mass of the monopole M (squares), and mass of the dual photon m 
(circles) vs. p. 
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An alternative method to demonstrate dual superconductivity is to detect London 
current in the flux tube configurations between qq 15 . For a detailed comparison of our 
approach with ref.16 we refer to ref.13. 

As for a direct determination of (11) it can be shown that (fx) has a gaussian 
distribution in the sense of the central limit theorem 13 . However, due to the exponential 
dependence on AS* the expectation value of (fx) is not centered at the minimum of AS. 
Let H{y) be the distribution of y = (3 AS — (3 (AS), the distribution probability for \x is 

PU) = II f ln(^) ] dln(^) 
V V J V 

If n(y) is gaussian with width a y , 

{n) = /Zexp((-^)) ct m = /Zexp((-^)) 

The fluctuations are larger than (//). 

A careful analysis requires the account for higher cumulants. The displacement of 
the maximum of P(p) with respect to (fx) should be kept into account when computing 
the so called constrained effective potential. 

As a final comment it can be shown that our construction gives the same result as 
that of ref.16 in the case of Villain action. 

We have thus a disorder parameter which is a reliable tool to detect dual super- 
conductivity. 

We have successfully repeated the construction for the X — Y model in 3d, where 
vortices condense to produce the phase transition. The result can be checked aganst 
experiment (liquid ife 4 ). 



3d X - Y model 21 (liquid #e 4 ) 

The model is defined on a 3d cubic lattice. An angle 0(i) is defined on each site i. 
The action reads 

S = ^^[l-co S (VW)] (54) 

i fi=0 

A^9(i) = 6(i + fx) — 9(i). The partition function is 

Z = JW^eMS) (55) 

i 

Z is a periodic functional of 9(%) with period 27r (compactness). Z is invariant under 
the change of variables 

- B{i) + f(i) (56) 

with arbitrary f(i) and so is any correlator of compact fields, in spite of the fact that 
the transformation (|56| ) is not a symmetry of the action. In running the indices in 
Eq.(|56|) from to 2 we anticipate that we shall consider the theory as the euclidean 
version of a 2+1 dimensional field theory. 
As (3 — > oo 

S ^ f£[A^)] 2 (57) 
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and the theory describes a massless scalar field. At (3 = /3 C = 0.454 a 2nd order phase 
transition takes place. Below (3 C vortices are expected to condense in the vacuum. 
Like for monopoles, condensation has always been demonstrated by the drop of the 
density of vortices when j3 raises through f3 c . We will show instead that a spontaneous 
symmetry breaking of the U(l) symmetry related to the conservation of vortex number 
takes place and that a legitimate disorder parameter can be defined. We define 

= <V (58) 

Under the transformation Eq. (^) A^ undergoes a gauge transformation 

—> + (59) 



The invariance under Eq. ([56]) means gauge invariance if the theory is phrased in terms 
of A^- From Eq.(|58|) it follows 

Fnu = dnA v - d v An = (60) 
In fact Eq.(|60D is valid apart from singularities. In terms of An 

Q(x) = f" expiiA^daf) (61) 
Jc 

and, if Eq.(|60D holds, the choice of the path C used in Eq.([H]) is irrelevant. A current 
in can be defined as the dual of 

J, = Sn a pd a A? (62) 

and 

d»Jn = (63) 

is the analog of Bianchi identities. The conserved quantity associated to Eq.(^) is the 
vorticity 

Q = [ d 2 xf(x, t)= f(VA A)q = j> Adx (64) 



/■ 

Since A = V8, single-valuedness of the action implies that Q = 2irn. If there are no 
singularities it follows from Eq.fl55|) that V A A = or, from Eq.(]64]) that Q = 0. 
There exist however configuration with non trivial vorticity. An example is 

0( q ) {x-y) = q arctan y^—^ (65) 

[x y)\ 

For this configuration 

A° = A=^-v e (66) 
r 

Vq being the unit vector in the direction (9) in polar coordinates. If A is the field of 
velocities 0( q ) describes a vortex. For this configuration j° = 2%8^(x — y) and 

A ■ dx = 27rq (67) 

c 

if the path C encloses y once. 

A disorder parameter describing condensation of vortices can be defined. In the 
continuum this is nothing but addition to the field 9(x) of any configuration of the field 
9 describing a vortex. The conjugate momentum to 9 being Psmd 9 = dC/d(d 9) 



fi(y, t) = exp 



i(3 / d x sm(d 9(x, t))9^(x — y) 



(68) 
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analogous to Eq. (|38|) . 

After Wick rotation the compactified version of Eq . (|68|) becomes (see Eq.(j37|)) 



p(y,t) = exp I -/3J2 [cos(A o 0(n, i) - 6 {q) {n - y)) - cos(A 6»(n, t) 



(69) 



When computing correlation functions of p(y, t) Eq.(|69|) produces a change of the term 
in the action containing A 8 at time t from A 9(n, t) to A 9(n,t) — 0( q ){n — y). The 
correlator 

x?(t) = (MM)Mo,o)> (ro) 

will behave as 

V(t) ~ (fj.) 2 + Ae- Mt (71) 

t — >oo 

and, from the definition of p 

Z[S + AS] 



V(t) 



Z[S] 



(72) 



The S + AS is obtained from S by the modification 

A 6(n, 0) -> A #(n, 0) - - y) (73) 

A #(n, t) -> A fl(n, t) + (n - $ (74) 

A change of variables 0(n, 1) — > 8(n, 1) — di q ) reabsorbs the modification in A o 0(n, 0), 
but shifts Ai8(n, 1) = Aj(n, 1) by a vortex Aj)^ = and sends 

A Q 0(n,l)^A o 0(n,l)-0( ? )(™-£> 

Similarly to what was done with £7(1) monopoles the construction can be repeated till 
n = t is reached and —9{ q ) cancels 9( g ) of Eq.fl731). 

T>(t) describes a vortex at y = propagating in time from to t. Instead of T>(t) 
p(t) = ^lnX>(£) can be studied. At large t 



pit) ~ p + Ae~ Mt 



(p) = exp 



1 

2 Jo 



Poo(P) d(5 



(75) 
(76) 
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From Eq.(|76"D the disorder parameter (p) can be determined. A typical behaviour 
of p is shown in fig. 6. In the thermodynamical limit (p) 7^ signals spontaneous 
symmetry breaking of the symmetry Eq. (|5"3"D , and hence condensation of vortices. At 
large f3 the theory is free and p can be computed by a gaussian integration. The result 
for a cubic lattice os size L is 



p = -11.33 • L + 72.669 (77) 

As L — > 00 this implies (p) = 0. For (3 < /3 C , p tends to a finite value, as L — > 00. 
Around /3 C a finite size scaling analysis can be done as in the U(l) model giving 

v = 0.669 ±0.065 [-670(7)] (78) 
0c = 0.4538 ± 0.0003 [.45419(2)] (79) 
5 = 0.740 ±0.029 (80) 



5 is the index of p. The numbers in parenthesis on the right are the accepted determi- 
nations by other methods 22 . The agreement is good. 

Monopoles in QCD. Revisiting the abelian projection 23 . 

At the classical level monopoles in non abelian gauge theories can be defined by 
the usual multipole expansion 24, 25 . At large distances the magnetic monopole field 
obeys abelian equations. By a suitable choice of gauge the direction of the Dirac 
string can be chosen, and magnetic charges are identified by a diagonal matrix in the 
fundamental representation with positive or negative integer eigenvalues. For SU(N), 
N — 1 magnetic charges exist, which correspond to a group {/(l)^ -1 . 

This classification coincides with the so called abelian projection. 

Let $<j be any field belonging to the adjoint representation: in what follows we 
shall refer to SU(2) for the sake of simplicity. For SU(N) the procedure is analogous 
with some formal complication. A gauge transformation U(x) which diagonalizes $<7 is 
called an abelian projection. U (x) will be singular at the locations where $ = 0. These 
locations are world lines of U(l) Dirac monopoles. U (1) is the residual invariance under 
rotation around the axis of $<? = $cr 3 after diagonalization. 

These monopoles are supposed to condense and produce dual superconductivity. 
There is a large arbitrariness in the choice of $<?, i.e. in the identification of monopoles. 
We will discuss this point in detail below. 

The meaning of the abelian projection can be better understood 26 in the Georgi 
Glashow model 27 . This is an 5*0(3) gauge theory coupled to a triplet Higgs 

C = -\g^G, v + (D^(D^) - V($) (81) 

with 

D f ,$=(d ll -gW fl A)$ 

and 

($ 2 ) 2 - /i 2 $ 2 " 

The model admit monopoles as soliton solutions, in the spontaneous broken phase 
where ($) = $° 7^ 0. 



V@) = 7 
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Usually a fixed (point independent) frame of reference is used in colour space, e.g. 
£j° (i = 1, 2, 3), the unit vectors of an orthogonal cartesian frame. A £° = e^kik- One 
can, however, define a Body Fixed Frame (BFF) by 3 orthogonal unit vectors £i(x) 

£i A £j = Sijk^k (82) 

such that £3(2;) = is parallel to the direction &(x) of the $ field. This system is 

defined up to a rotation around £3. An element of the gauge group exists R(x) such 
that 

Ux) = R-\x)S (83) 
By construction R(x) is the gauge transformation which operates the abelian projection, 
since R(x)£ 3 = £3. Since £^ = 1 

<9 M & = ^ M A ii 

or 

D„(u;)£ = (0„ - u^A)ii = (84) 

Eq.(|84|) implies [D^, D u ] = G^iQ) = 0. The field oj^ is a pure gauge. The last equality 
reads explicitely 

dpUv - dvUjp + A uj u = (85) 
Eq . (|85|) is true apart from singularities. As a consequence of Eq.(|8~5"D 



$(x) = Pexp 



Q.-Tdx 11 

c 



$° (86) 



Due to Eq. (|85| ) the integral in Eq. (^) is independent of the path C. This is true apart 
from singularities which can give a non trivial connection to space time. 

A t'Hooft 28 Polyakov 29 monopole configuration has a zero of $(a;) at the location 
of the monopole, and in that point the abelian projection R(x) has a singularity. 

The singularities of cJ M can be studied by expressing £(x) in terms of polar coor- 
dinates 6(x), ip(x), with polar axis 3 in colour space. The singularities come from the 
fact that i[)(x) is not defined at the sites where 9(x) = 0,7r. In terms of G^ v {oj) the 
potentially singular term at 9{x) = 0, tc is 26 

F 3 , u (x) = -cosd(x)(d,d u - d u d^(x) (87) 

The singularity exists where 9 = 0, 7r or at the sites where $ is in the direction 3, and 
the field is parallel to $ in colour space and abelian. The singularity is a string with 
flux ±2n7c (n = 0,1...). The abelian field is the field related to the residual U(l) 
symmetry along the 3d axis. The field strength is the abelian part of F^, or 

J> = $ ■ G,u ~ • {Dfi A (88) 
9 

Indeed, in the abelian projected frame = gA^ A £3 and 

-$(D M $ A D v $) = g{W^ A W u )$ 
9 



which cancels the non abelian term of ^G^ in Eq. (p8|) . is nothing but a covariant 
expression for the U(l) field identified by the abelian projection. 

is a gauge invariant quantity. For a t'Hooft Polyakov monopole configuration 

is the field of a pointlike Dirac monopole, located at the zero of $(x). 

In QCD there are no fundamental Higgses. The idea is that monopoles could be 
defined by any composite field $(x) in the adjoint representation 23 . No unique criterion 
is known for the choice of the operator $ in QCD. Popular choices are 23 
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1) is the Polyakov line 

r x o -, -, 

P(x,x )=f P exp(L4 M • T dx^) (89) 

J C,xo 

The path C being the line x = constant along the time axis closing at infinity by 
periodic boundary conditions. 

2) Any component F^ u of the field strength tensor. 

3) The operator implicitely defined by the maximization of 

£Tr{^>y 3 f/>)} (90) 

This choice is known as maximal abelian projection. 

4) For SU(3) F 2 U since, contrary to the case of SU(2), it is not a singlet, due to the 
d algebra. 

For each of these choices a U(l) gauge field (actually 2 for SU(3)) can be identified, 
which couples to monopoles, and a creation operator of monopoles can be constructed 
on the same lines as for compact U(l). 

A strategy to answer the question whether QCD vacuum is a dual superconductor 
is to detect the condensation of the monopoles defined by different abelian projections 
in the confined phase and across the transition to the deconfined phase. As shown in 
the previous sections a reliable tool (disorder parameter) exists for that, which has been 
successfully tested in systems which are well understood. Before proceeding to that we 
will test the ideas of this section on a simple system: the Heisenberg ferromagnet. 

The Heisenberg ferromagnet 30 . 

The action is 

£ = ^EE \n(x)A,n(x) = £ £(1 - *(* + ( 91 ) 

x runs on a cubic 3d lattice and n is a vector of unit length, n 2 — 1. 

We shall look at the model as a 2+1 dimensional field theory. The Feynman 
functional 

Z = I Y[<m(x) exp(-pS) (92) 

J X 

has a much bigger symmetry than the lagrangean. Any local rotation n(x) — > R(x)n(x) 
even if it does not leave C invariant, is reabsorbed in a change of variables in the Feyn- 
man integral, leaving the measure invariant. Assuming constant boundary condition 
at infinity, n , we can write 

n{x) = R~\x)n G (93) 

R^ 1 is determined up to a rotation along n. As in the Georgi Glashow model a gauge 
field cJ M can be defined by introducing a Body Fixed Frame with £ 3 = n. Then 

d^i = ijJfx A £i 

or 

D^i = {dp - \u>fi% = 
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which implies F^ v {uo) = 0, apart from singularities and 



n(x) = Pexp 



/ 

Joo,C 



■Tdx^ 



n 



(94) 



independent of the path C, being a pure gauge. Again this is true apart from 
singularities. A conserved current exists 



07T 



(95) 



jfi is parallel to n, since both d a n and dpn are orthogonal to it. The corresponding 
conserved quantity is 



Q = J d 2 xj°(x) = ^J d 2 xn ■ (d 1 nAd 2 n) 



(96) 



which is nothing but the topological charge of the 2 dimensional version of the model. 
Instantons of the 2 dimensional model look as solitons of the 2+1 dimensional one. 
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— * 

can be also written as 



n 



07T 



(97) 



Except for the singularities corresponding to the locations of the instantons £u 1 A uj 2 = 
d\UJ2 — d 2 oJ\ and 

Q = —j>Q.fidx (98) 

A non trivial connection is created by the presence of solitons. A direct calculation 
shows that the field has Dirac string singularities propagating in time from the 
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center of the soliton. The transition from a magnetized phase to a disordered phase 
can be seen as a condensation of solitons. 

The system undergoes a second order phase transition at f3 = f3 c ~ 0.7 from the 
magnetized phase to a disordered phase. We have investigated if the solitons described 
above condense in the disordered phase. A disorder parameter can be constructed, on 
the same line as in the compact U(l) and in 3d X — Y model, as the vev of the creation 
operator of a soliton, as follows 

where S + AS is obtained from S by the change 

(A n(x , x )) 2 = (n(x, x + 1) - n(x, x )) 2 -> (R^n^x, x + 1) - n(f , x )) 2 (100) 

R q is a time independent transformation which adds a soliton of charge q to a configu- 
ration. 

Numerical simulations show that in the thermodinamical limit (p(x)) vanishes in 
the ordered (magnetized) phase, is different from zero in the disordered phase, and at 
the phase transition obeys a finite size scaling law from which the critical indices and 
the transition temperature can be extracted. A typical form of p = d\n(p) / d(3 is shown 
in fig. 8. 

The results are still preliminary but agree with the values known from other 
methods 31 . We get 

P c = 0.69 ± 0.01 (101) 

and 

z/ = 0.7±0.2 (102) 

This shows that the phase transition to disorder in the Heisenberg magnet can be 
viewed as condensation of solitons. 

The string structure of the singularities of the field cJ M is similar to Georgi Glashow 
model, and the field strength tensor is generated by the topology of the field n, even in 
the absence of gauge fields. Indeed going back to the previous section, the monopoles 
only depend on the Higgs field. Monopoles exposed by the abelian projection are not 
lattice artifacts, but reflect the dynamics of the field $. 



Dual superconductivity in QCD. 

I will present the first results of a systematic exploration of monopole condensation 
below and above the deconfining transition in different abelian projections 33 ' 32 . 

Besides the disorder parameter, we also measure, at T = 0, the monopole anti- 
monopole correlation function to extract a lower limit to the effective Higgs mass, as 
well as the penetration depth of the electric field, i.e. the mass of the photon which 
produces (dual) Meissner effect. We use for that APE QUADRIX machines. 

The creation operator for a monopole is constructed in analogy with the U(l) 
operator as follows: the II * plaquettes in the action at the time, say n = when the 
monopole is created are modified as follows 

U 0i (n,n ) = Tr{Ui(n,0)U o (n + i,0)U^(fi,l)Ui(n,0)} 

- Tr {ul(n,0)U o (fi + i,0)U}(fi, l)uUn,0)} (103) 
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In Eq.(|T03D 

U'(n, 0) = e iA ^ exp(i^a^^)U t (n, 0) exp(i<l>o ; ^|^)e- iA(n+1) 

We call S + AS the resulting action. Then 

Z[S + AS] 



00 



Z[5] 



(104) 



(105) 



The vector potential describing the field of the monopole 6, has been split in a transverse 
part bf- with <9jfej = and a gauge <9jA. The gauge dependence in Eq.( [L05| ) can be 
reabsorbed in a redefinition of the temporal links, which leaves the measure invariant. 

A redefinition of C/j to £7 = e 1<73 2 [7e 1<T3 2 in the abelian projected gauge can 
be reabsorbed in a change of variables which leaves the measure invariant. The space 
plaquettes however acquire in each link factors 



&■ (n) _ « b^(n) 



In the abelian projected gauge the generic C/j canbe written as 



(106) 



(107) 



#i = a, + exp(io"3^j) is the abelian link. The transformation Eq.( |104|) adds hi to $i 
or adds the magnetic field of the monopole to the abelian magnetic field 



AiOj - AjOi -> Ai9j - AjOi + Aibj - Ajb, 



(108) 



An abelian projected monopole has been created, (/jl) or better p = d{hx{p)) / dp is 
measured across the deconfining phase transition. Typical behaviours are shown in 
fig.9 and fig. 10 for the monopole defined by the Polyakov loop. 
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Fig.9 p v.s. P 577(2) gauge theory. (Lattice 
12 3 x 4) 
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A careful analysis of the thermodynamical limit produces evidence of dual supercon- 
ductivity below the deconfining temperature T c , for different choices of $, Polyakov line, 
F^v, max abelian. 

We need more work and more statistics to determine the critical indices and the 
type of superconductor. 



Discussion 

Most of the success of the dual superconductivity mechanism for confinement is 
presently based on the abelian dominance and on the monopole dominance numerically 
observed in the maximal abelian projection 34, 35, 36 . In SU(2) gauge theory after max 
abelian projection, which amounts to maximize the quantity 

A = 52Tr{a 3 UM<T 3 Ul(n)} 

with respect to gauge transformation, all the links are on the average parallel to 3 axis 
within 80 — 90%. The string tension computed from the abelian part of the Wilson loops 
accounts for 80 — 90% of the full string tension. Similar behaviour is found for many 
other quantities. This is called abelian dominance. In addition the abelian plaquettes 
can be split as in Eq.(|29|) in a monopole part, n^ u and the residual angle 8^, usually 
called Coulomb field. The empirical observation is again that the contribution of the 
Coulomb part to the abelian quantities is a small factor, so that, for example, the string 
tension is dominated in fact by the contribution of the abelian monopoles (monopole 
dominance). 

Apparently such dominance is not observed in other abelian projections, except 
maybe in the Polyakov line projection, where, if the string tension is measured from 
the correlation between Polyakov lines it is 100% abelian dominated by construction. 

The more or less explicitely expressed idea is then that some abelian projection 
(the max abelian) is better than others and is ' 'the abelian projection", identifying 
the degrees of freedom relevant for confinement. On the one hand the fact that in 
this projection links are in the abelian direction at 80 — 90% implies dominance as a 
kinematical fact: on the other hand maybe it is non trivial that such a gauge exists. 
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The attitude presented in this paper is somewhat complementary: dual supercon- 
ductivity is related to symmetry, and the way to detect it is to look for symmetry. From 
this point of view, independent of possible abelian dominances, what we find is more 
similar to the idea of t'Hooft that all abelian projections are physically equivalent. We 
find condensation of max abelian, Polyakov line, field strength monopoles. 

There are few aspects of the mechanism which need further understanding. 

1) Any abelian projection implies that the gluons corresponding to the residual 
£7(1) 's have no electric charge with respect to them, and hence are not confined. 
This contrasts with the observation reported in sect. 1, about the string tension 
in the adjoint representation. 

2) If the mechanism of confinement were superconductivity produced by condensa- 
tion of the monopoles defined by some abelian projection then the electromagnetic 
field in the flux tubes joining q — q pairs should belong to the projected U(l). 
Lattice exploration show that it is isotropically distributed in color space 6 ' 37 . 

3) There are infinitely many abelian projections which can be obtained from each 
other by continuous transformations, e.g. by shifting the zero of $>(x). If one of 
them were privileged, it is hard to understand how others, which differs by small 
continuous changes, could be so different. 

Most problably a more complicated and new mechanism is at work, a kind of non 
abelian dual superconductivity, which manifests itself as abelian superconductivity in 
different abelian projected gauges. Our aim is to understand it. 

It is a special pleasure to thank all the collaborators who contributed to this 
research, in particular Giampiero Paffuti, Manu Mathur and Luigi Del Debbio. 
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